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Abstract 

The search for heavy resonances in the dijet channel is part of the on-going physics 
programme, both at the Tevatron and at the LHC. Lower limits have been placed 
on the masses of dijet resonances predicted in a wide variety of models. However, 
across experiments, the search strategy assumes that the effect of the new particles 
is well- approximated by on-shell production and subsequent decay into a pair of jets. 
We examine the impact of off-shell effects on such searches, particularly for strongly 
interacting resonances. 
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1 Introduction 

Dijet production is an integral part of studies conducted at any hadron collider. Apart 
from being important for the ratification of our understanding of QCD, this process also 
provides a fertile ground for new physics searches, particularly searches for new particles 
that may appear as resonances in the dijet invariant mass spectrum. The large production 
cross-sections mean that significant conclusions may be drawn even with low integrated 
luminosities. Of course, QCD itself gives rise to a large background which is a challenge that 
any new physics search has to contend with. However, in case the new physics signal appears 
in the form of a localised excess in production rates in a particular region of phase space 
{e.g. a resonance), it is often possible to subdue the background with appropriate kinematic 
cuts. 

While no new physics beyond the Standard Model (SM) has been discovered yet, exper- 
imental collaborations at both the Tevatron and the LHC have been using the dijet channel 
to place limits on masses of new particles in a wide variety of models [H El [3] . The strategy 
has been to compute upper bounds on a ■ B ■ A (where a is the cross-section for on-shell 
production of the particle, B is the branching fraction into a pair of jets and A is the ex- 
perimental acceptance) and hence rule out certain mass ranges for the new particles. In 
essence, a ■ B is the excess that arises due to the presence of the new particle. This implies 
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the assumption that the new particle under consideration has a small width. However, for 
some strongly decaying particles, which the Tevatron/LHC have excluded upto a large mass, 
this assumption may not hold true. Examples are particles such as axigluons and colorons 
which typically have a width > 10% of their mass and, are certainly not 'narrow resonances'. 
However, in the context of a hadron collider, it is also important to remember that the 
'observed width' as would be measured from the invariant mass distribution is not the true 
decay width of the particle associated with the resonance. The shape of the resonance gets 
distorted by fluxes of the initial state particles which, in turn, depend on the parton level 
center-of-mass energy. Further, experimental reconstruction of jets is a challenging task. 
Limitations of detector resolution and reconstruction algorithms affect the resolution of the 
jet-jet invariant mass rrijj. If the aforementioned effects overwhelmingly dominate the ob- 
served width, then the theoretical approximation of 'narrow-width' may be inconsequential. 
In this note, we examine such aspects of dijet resonance searches in an attempt to compare 
the relative merits of different search strategies. 

The rest of this article is organised as follows. We begin by briefly recapitulating, in 
Section a few model templates with particular emphasis on their status vis a vis searches 
at the Tevatron and the LHC. Section [3] details our calculations and delineates the numerical 
effects due to the off-shell contributions. Finally, we conclude in Section |H 

2 Dijets and new physics 

With both the Tevatron and the LHC accumulating substantial luminosity, several stud- 
ies have looked at possible resonances. While resonances decaying into leptons (such as a 
hypothetical Z') are relatively easy to look for, of particular interest are particles that prefer- 
entially decay into hadronic states. Indeed, considerable interest has been generated by two 
Tevatron measurements, viz. the reported excess [1] in the Wjj final state on the one hand, 
and the longstanding discrepancy in the forward-backward asymmetry in tt production [5] 
on the other. While several models have been proposed as solutions to these deviations from 
the SM, it is imperative that they be examined vis d vis other processes. With each such 
explanation positing new couplings involving quarks, dijet production is a natural theatre 
for such investigations [Hj. Indeed, a recent study [7] has attempted a model independent 
study of colored resonances. 

The simplest algorithms for new particle searches naturally concentrate on situations 
wherein the role of the said particle, in any process involving only SM particles as asymp- 
totic states, can be well-approximated by a narrow resonance. Unfortunately though, this 
approximation is often rendered invalid. Further, under certain conditions, even the assump- 
tion of a constant width may need correction. Numerous examples abound in low-energy 
hadron physics (see for example, Ref. [8]). Indeed, this effect turned out to be significant even 
for a narrow resonance such as the Z [5]. More recently, this question has been examined in 
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the context of production of top-pairs [TOl [TT], dijets [12], a heavy SM Higgs boson [13] and 
hypothetical W's [H] and Z'sQ [B]. 

In this study, we consider axigluon^ and colorons as templates for broad dijet resonances. 
These are color octet gauge bosons appearing in certain classes of models that hypothesize 
an extended color gauge group SU{?>)a ® SU{3)b at high energies. The extended gauge 
symmetry is broken spontaneously to SU{3)c that one associates with strong interactions 
in the SM. Axigluons and colorons correspond to the broken generators in the respective 
models and, hence, are massive. The coupling of axigluons to quarks is given by gs1ij.'j5 
while that of colorons is (7^7^ cot ^, ^ being the angle that characterizes the mixing between 
the two SU{3) groups. In certain variants of the model, the coloron couples preferentially to 
the third generation quarks. However, here we consider only the flavor-universal coloron, for 
which cot ,^ < 4 in order that the model remains in its Higgs phase. Details of axigluon and 
coloron models can be found in Refs. [18], [19] and Refs. [201 EI] respectively. Over the years, 
limits on axigluon and coloron mass (denoted henceforth by Ma and Mc, respectively) have 
been upgraded continuously based on experiments that have been in operation at various 
times [22]. The current experimental limits are due to searches in the dijet channel at the 
LHC. The ATLAS experiment has ruled out Ma< 3.32 TeV [5] while CMS has ruled out 
axigluons and colorons of mass below 2.47 TeV [2]. Before this, the CDF experiment at the 
Tevatron had placed a lower limit of 1250 GeV on Ma, Mc [I]. However, barring the analysis 
of Ref. [To] (included by the Particle Data Group [23]), the rest implicitly incorporate the 
narrow-width approximation, and thus need to be accepted with care. 



3 Numerical Analysis 

In this analysis, for the Tevatron as well as the LHC, we consider representative masses of 
the aforementioned particles in the range that is well within the kinematical reach of the 
respective machines. Then, for each such representative case, we compute the signal in two 
ways - first, by considering on-shell production and subsequent decay of the new particle, 
and second, by computing the full cross-section including off-shell effects and all contributing 
amplitudes (s-, t- and u- channel etc.). We label these two cases by OS and FL respectively. 

The computations are performed using CalcHEP [21]. The CTEQ6L parton distribu- 
tions [25] are used along with the compatible value for as as obtained using the 2-loop 
/S-functioiJl. The renormalization and factorization scales are set to the sub-process center 
of mass energy. 

^Although particles such as a heavy SM Higgs, W or Z' do not decay strongly, their large masses 
(> 600 GeV) still result in a substantially large width. Of particular importance is the fact that, for such 
widths, the interference with non-rcsonant contributions to the amplitude become important, thereby ren- 
dering the on-shcU approximation rather inaccurate. 

■^Axigluons arc of particular interest in the context of the forward-backward asymmetry in production. 
See, for example, Refs. [I0l[li[l7]. 

^While this might seem paradoxical given that we are computing only the leading-order matrix elements, 
note that the CTEQ collaboration uses as calculated at NLO to extract the CTEQ6L distributions. Further- 
more, the use of CTEQ6L1 distributions alongwith as(LO) does not qualitatively change our conclusions. 
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3.1 At the Tevatron 



As a representative example close to the lower bound achieved by CDF, we consider an 
axigluon of mass Ma= 1.2 TeV. The natural width of this particle is 100 GeV. To enhance 
the signal to background (SM) ratio, we only generate events where the jets have a pj- of at 
leasto 150 GeV. Since we are primarily interested in a localised excess in the rrijj distribution, 
we must effect a fit to the rrijj spectrum significantly away from the excess. An excellent fit 
is obtained [1] in terms of the four-parameter function f{x) = ao (1 — x)""^ a;a2+a3 inx -y^^j^gj-g 
X = rajjl ^fs^. Interestingly, a Gaussian in x gives a (three-parameter) fit that is virtually 
as good (in terms of per degree of freedom). 

Once the aforementioned (theoretical) SM spectrum is obtained, the invariant mass spec- 
trum for both the OS and FL case can, then, be compared with it. In the infinite resolution 
limit, the OS distribution would be characterised by a sharp spike in the bin corresponding 
to rrtjj = 1.2 TeV. However, in an experimental situation, the spike gets smeared due to 
detector resolution effects. Bearing this in mind, we apply a Gaussian smearing to the en- 



ergy of all the final state particles with 5Et/Et = 50%/y -E'T(GeV) © 3% (this being the 
resolution of the central hadron calorimeter for CDF [1]). To account for possible upward 
scaling of the energy, we now impose pt > 200 GeV on the two leading jets. The resultant 



rrijj distribution is plotted in Fig, 1(a) The difference between the OS and the FL cases is 
clearly visible even to the naked eye and this shows that, even after smearing, the difference 
between the two distributions is discernible. 



It could be argued that, in making Fig, 1(a), we have taken the liberty of finely binning 



the events. While this demonstrates the difference between the two situations very well, such 
small bin sizes may not be achievable in practice. The study of Ref. [1], for example, uses 
a bin width of O.lrrijj. To examine the effect of this reduced resolution, we redistribute the 



events accordingly. The result is shown in Fig, 1(b) 



Understandably, the difference between the FL and the OS approximation is reduced to 
an extent. Nonetheless, a significant difference between the two does persist as it does with 
the SM. At this point, one must ask how significant the difference is statistically. To this 
end, we calculate the binwise significance Si defined as 

<Si = Si/ \f^i , 

where Si [Bi) are the number of signal (background) events in a given The results, 



plotted in Fig, 2(a), bring out two facts. First, near the resonance, the significance obtained 
is much larger for the OS case. Larger expectations of significance imply greater sensitivity 
and allow the model to be ruled out at higher confidence levels. However, in this case, 
the greater significance is but an artefact of the narrow-width approximation which, as 
demonstrated by Fig ]l(a)| is clearly not applicable in this case. Second, away from the 
resonance, cross-sections are suppressed and are /ower than the SM prediction if one considers 
the full (FL) amplitud^. This is an effect of the interference of the multiple axigluon- 



^This is nothing but the imposition of CKIn(3) > 150 GeV in Pythia [26]. 

^With the rebinning, the number of events in each bin is large enough for the Poisson-distributed number 
of events to be well-approximated by Gaussians. 

^This suppression was not visible in Fig 1(a) simply on account of the scale of the plots. 
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Figure 1: The dijet invariant mass spectrum at the Tevatron corresponding to Ma= 1200 GeV 
(the corresponding width is 100 GeV). An integrated luminosity of 10 fb~^ has been 
assumed and a restriction of px > 200 GeV imposed on each of the two highest p^ jets. 
The left panel corresponds to a constant rrijj bin width, while the right panel redistributes the 
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mediated contributions (including t- and u-channel diagrams) with the SM amphtudes. In 
fact, with py*" = 200 GeV, the total cross-section is slightly lower than the corresponding 
SM prediction. While the dominant new physics contributions pertain to qtqi — )■ g^gj, even 
subdominant processes such as qiqj — )■ qiqj are suppressed on account of the destructive 
interferences engendered by the t-channel axigluon exchange contributions. Understandably, 
no such suppression exists in the OS approximation. 

It is interesting to contemplate the impact of kinematic cuts on this difference. It seems 
plausible that placing a strong cut on the jet rapidities would not only eliminate a larger 
fraction of the SM background, but also significantly reduce the contribution from the in- 
terference of the axigluon amplitudes with the t-channel gluon-exchange ones. In effect, the 
imposition of stronger requirements on the jet transverse momenta would be expected to 



push the new physics contribution closer to the OS approximation. In Fig, 2(b) we show the 



results for the significance on demanding = 500 GeV and \y\ < 0.7. Comparing it to 



Fig, 2 (a) , the increase in significance levels is obvious. What is even more striking is that the 
increase in significance is more pronounced for the OS approximation than it is for the full 
(FL) calculation. In other words, even the imposition of such strong cuts does not validate 
the OS approximation, and the large sensitivity is partly illusory. The absence of a dip in the 



differential cross section (as compared to Fig 2(a)) is understandable as the corresponding 
rrijj bins have been eliminated by the imposition of the strong px cut. 

Further, if one considers the background subtracted excess (as shown in Fig. [3]), one finds 
that, in case of OS, the excess may be approximated by a Gaussian distribution. On the 
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Figure 2: Binwise significance expected, at the Tevatron, for an axigluon of mass 1.2 TeV 
and bin-widths ofnijj/lO. The two panels correspond to differing requirements on the two 
leading jets (a) pr > 200 GeV, and (b) pr > 500 GeV, \y\ < 0.7. 



other hand, the FL distribution is asymmetric and is not amenable to fitting by a simple 
functional form. Quite naturally then, if one were to attempt an extraction of parameters 
such as the mass and width, FL and OS would give rise to different values. In fact, this is 
borne out quite clearly by Fig. [31 Not only is the position of the peak visibly shifted in case 
of FL, the shape too is sufficiently skewed to disallow an agreement with an OS Monte Carlo. 

While the above discussion involved the invariant mass distribution which is the crucial 
observable in any resonance search, it is possible to appreciate the difference between OS 
and FL even by considering just the deviation in the total cross-sections, albeit with strong 
kinematic cuts. In FigJUwe show the deviation from the SM dijet cross-section as a function 
of Ma. It is clear that considering simply on-shell production results in over-estimation of 
the signal for most cases that are within the kinematic reach of the Tevatron. Again, this 
is a consequence of the multiple interference terms involving gluon/axigluon mediated s-, t- 
and M-channel amplitudes. In fact, were one to relax the cuts (for example, imposing only 
p-mm _ GeV with no restriction on rapidity), the effect of the interference would drive 
the total cross-section below the SM cross-section. It should be noted though that claiming 
a discovery (or otherwise) based on total rates alone is fraught with danger as this observable 
is particularly sensitive to higher order corrections, uncertainty in parton distributions etc. 

The case of the coloron is even more curious. The coupling strength is Qs cot C,- With 
1 < cot,^ < 4, the widths could easily be much larger than that for an axigluon of similar 
mass (~ 30% of the mass with cot,^ of just 2). As a result, it would be very difficult, 
if not impossible, to identify the excess as a resonance (see FigJS]). Although it might be 



6 




1050 1100 1150 1200 1250 1300 1350 1400 
mjj (GeV) 

Figure 3: The excess at the Tevatron (assuming Ma= 1200 GeV and with = 200 GeV) 
after subtracting the SM background for the two scenarios FL and OS. The curves depict the 
result of the fitting when a Gaussian distribution is assumed. 
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Figure 4: The excess dijet cross section at the Tevatron as a function of the axigluon mass 
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argued that the enhancement in rates in the high rrijj regioiu, in itself, would constitute a 
smoking gun signal, note that the task is not as straightforward. The wide region (in the 
rrijj spectrum) of the coloron's influence means that virtually no part of the spectrum can 
be termed to be essentially SM-like, thereby changing the entire nature of the fit algorithm. 
Thus, a much more sophisticated algorithm, including a higher order calculation of the SM 
dijet spectrum, would be necessitated; in particular, such corrections can, and do, change 
the shape of the spectrum. 
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Figure 5: The dijet invariant mass distribution, at the Tevatron, in the presence of a coloron 
with Mc= 1200 GeV and cot^ = 2 (this translates to Tc= 408.5 GeV). Each of the two 
leading jets is required to have pr > 200 GeV and an integrated luminosity of 10 fb~^ has 
been assumed. 



3.2 At the LHC 

At the LHC, the pitch is queered by the low fluxes for antiquarks in the initial state. Due to 
this, the contribution from the s-channel amplitude suffers a reduction as compared to the 
contribution from the t- and w-channel amplitudes which may have two quarks in the initial 
state. Consequently, the OS approximation is less likely be valid here as compared to the 
case of the Tevatron. For the SM background, a fit analogous to that for the Tevatron can 
be made [2], of course with differing parameters. Once again, a Gaussian fit works almost 
as well. 

^In fact, one of the motivations behind the proposal of the flavor-universal coloron model was the excess 
reported by the CDF collaboration in the tail of the i^T-spectrum for inclusive jet production [57]. However, 
when more data was analysed, this discrepancy disappeared 
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With the mass reach of the LHC being greater, we consider Ma= 2.0 TeV (rA= 160.2 
GeV), and, to enhance the signal to background ratio, we impose cuts much stricter than 
those for the Tevatron, namely = 700 GeV and \y\ < 0.7. For the energy resolution, we 
assume 6Et/Et = 70%/ ^'^(GeV) © 8% [29]. After smearing, a resonance-like structure 



seems apparent for the OS case, though no such claim can be made for FL (Fig. 6(a)). 
However, once larger (and experimentally more reasonable) bin-width^ are considered, the 
difference is expected to reduce. Note that part of this effect is offset by the increased 



statistics in each bin. As Fig. 6(b) shows, even with larger bin-widths, OS yields higher 
values for significance close to the resonance, whereas FL gives rise to significant deviations 
away from the resonance. 
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Figure 6: (a) T/ie dijet invariant mass spectrum at the LHC corresponding to Ma= 2.0 TeV 
(rA= 160.2 GeV). An integrated luminosity of 10 fb~^ has been assumed and a restriction 
of pt > 700 GeV and \y\ < 0.7 has been imposed on each of the two highest pt jets, (b) 
Expected binwise significance for the same case. 

Finally, we present the comparison between the invariant mass distributions for a coloron 
with Mc= 2.0 TeV and cot^ = 2 (rc= 645.7 GeV) in FigJH As with the case of the 
Tevatron, the OS approximation is demonstrably a very poor one and almost the entire dijet 
mass spectrum gets modified. 



3.2.1 Other Observables 

We have, until now, concentrated only on the invariant mass distribution as a discriminator, 
not only between the SM and new physics, but also amongst various forms of the latter. With 
the resolving power decreasing as the natural width increases, the case for other variables 
such as angular distributions becomes progressively stronger. For example, it has been 

^Though we use bin-widths of O.lnijj here, at ATLAS, for example, a resolution of O.OSmjj is possible [50] . 
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Figure 7: The dijet invariant mass spectrum at the LHC corresponding to Mc= 2.0 TeV and 
cot^ = 2 (rc= 645.7 GeV). An integrated luminosity of 10 fb~^ has been assumed and a 
restriction of pt > 700 GeV and \y\ < 0.7 has been imposed on each of the two highest px 
jets. 



argued [12] that the latter could be used to distinguish between a broad resonance and 
contact interactions. 

In this context, we re-examine a variable that has been used by the ATLAS collabo- 
ration [31] to obtain limits on quark contact interactions. In essence, they compare the 
number of events where both jets satisfy \y*\ < 0.6 (where y* is the rapidity in the partonic 
CM frame) with that where the extent of centrality is relaxed to \y*\ < 1.7. To be precise, 
F-^ijrijj) is defined as [31] 

, . _ Nl,,^,M\ < 0-6) 

where, the superscript i denotes the i^^ bin in the mjj distribution. We compute the 
distributions for an axigluon of mass 2000 GeV for each of the FL and OS cases. As suggested 
by Ref . [31] , the distribution is indeed sensitive to mass dependent changes in production 
rates in the central rapidity region. But that is not all. It is also sensitive to the assumption 



of a narrow- width approximation. In other words, as can be seen in Fig. 8(a) the shape of 
the distribution is markedly different in the two case^. The difference persists even when 



the bin-size is increased (Fig, 8(b)) 



^That the shape of our SM curve are different from that in Ref. [31] is attributable to the strong pT cut 
imposed by us. 
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Figure 8: The F^{mjj) distribution at the LHC for Ma= 2.0 TeV with p^™= 500 GeV. In 
(b), the bin-width taken to be O.lnijj. 

4 Summary 



New physics scenarios often predict resonances in the dijet channeL Very often, search strate- 
gies focus on setting hmits on a ■ B, the product of the cross section for on-shell production 
and the branching fraction into the dijet channel. While this strategy is perfectly valid in 
many cases, its applicability is questionable in the case of particles which have large widths 
due to either strong interactions or large masses. It has been argued in the literature [1], 
though, that the natural width plays only a subservient role to detector resolution effects, 
and, consequently, the former is oft neglected in experimental strategies. 

To examine this assertion in a qualitative manner, we considered two models and, for each, 
performed a simple analysis comparing the signal profile with and without the 'narrow-width 
approximation'. The dominant experimental effects (resolution) were incorporated through 
a smearing of the energies of the final-state particle^lj. We do find that the profiles can be 
substantially different in the two cases. For example, in the axigluon case, the natural width 
broadens the invariant mass distribution to a significant extent over and above the broadening 
due to resolution. A further complication arises from the fact that, in addition to the s- 
channel diagrams, we also have t- and w-channel contributions. While these can be safely 
neglected for a narrow resonance, their importance increases as the width becomes larger. 
In addition, channels such as ud — )■ ud which do not proceed through s-channel diagrams 
do receive contributions from a t-channel axigluon exchange. In totality, it transpires that 
inclusion of such effects would lead to tangible worsening of the exclusion limits. 

^"Indeed, we made a conservative choice of resolution parameters so as to enhance the experimental effects. 
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At the LHC, these effects could be even more important. The reason is not difficuh to 
fathom. At the LHC, the g is a sea-quark and, consequently, its density is relatively small 
for high Bjorken-a;, the regime that producing a very heavy resonance would require. On 
the other hand, t-channel axigluon contributions to qq — > qq remain unsuppressed. Thus, a 
careful analysis is even more mandatory in that arena. 

The signal to background ratio is, of course, dependent on the kinematical cuts. We have 
demonstrated that this dependance itself is quite different for the OS case as compared to the 
full calculation. Thus, we feel that the analyses should preferably be done with the inclusion 
of all diagrams and, possibly retuning the selection cuts so as to obtain robust exclusion 
limits. Moreover, while we have presented only tree-level results here, contributions from 
higher orders in perturbation theory would also have to be taken into account, both for the 
SM as well as the relevant new physics scenario. 

While the preceding discussion has largely focussed on the axigluon, it is by no means 
the only example of a colored resonance with a large width. Indeed, colorons generically 
have even larger widths. As we have shown explicitly, considering even cot^ = 2 (where, 
theoretically, 1 < cot ^ ^ 4) would render the width quite large and would almost totally 
obliterate a resonance. Rather, a very broad excess would be visible. Such an excess, in 
principle, could stand out once the background has been subtracted. However, the very 
width makes it virtually impossible to fit a functional form to the observations away from 
an excess, for, indeed, the excess spans almost the entire available invariant mass range. 
Hence, identification of the excess would require either a very precise knowledge of the SM 
prediction and perhaps even new experimental techniques. 
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